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ABSTRACT 

Superluminal motion is a common feature of radio jets in powerful 7-ray emitting active galactic nuclei. 
Conventionally, the variable emission is assumed to originate near the central supermassive black-hole where 
the jet is launched on parsec scales or smaller. Here, we report the discovery of superluminal radio features 
within a distinct flaring X-ray emitting region in the jet of the nearby radio galaxy M87 with the Very Long 
Baseline Array. This shows that these two phenomenological hallmarks - superluminal motion and high- 
energy variability - are associated, and we place this activity much further (> 120 pc) from the "central engine" 
in M87 than previously thought in relativistic jet sources. We argue that the recent excess very high-energy TeV 
emission from M87 reported by the H.E.S.S. experiment originates from this variable superluminal structure, 
thus providing crucial insight into the production region of 7-ray emission in more distant blazars. 
Subject headings: Galaxies: active — galaxies: jets — galaxies: individual (M87) — radio continuum: galaxies 
— radiation mechanisms: nonthermal 



1. INTRODUCTION 



The proximity of M87 (Z)=16 Mpc; lTonrvlll99lb makes it 
one of the best systems to study relativistic jets at high linear 
resolution (Fig. [T]). Our observations with the Chandra X- 
ray Observatory is olated short (month) timescale variability 
dHarris et al.ll2006h in ajet region previously dubbed 'HST-l', 
which is separated by 0.86" (60 pc, projected) from the cen- 
tral "active" supermassive black hole (SMBH). The variability 
culminated in a factor >50 X-ray outburst making HST-l the 
brightest X-ray source in the galaxy for a few years. Figure|2] 
shows a lightcurve from 2000 to the end of 2006. 

Observations with the Hubble Space Telescope (HST) and 
the Very Large Array (VLA) show comparable activity in 
the optical and a t radio frequencies (Perlman et alJ 120031 ; 
[Harris et"al]l2006h . After the detection of appreciable radio 
flux in HST-l from our first season of VLA observations (in 
2003), we began to monitor the jet at higher (sub-parsec) res- 
olution with the NRAO^ Very Long Baseline Array (VLBA) at 
three frequencies (0.33, 0.61, and 1.7 GHz) commencing in 
Jan. 2005. This was the period at which the X-ray and opti- 
cal intensities were peaking, and the radio intensity plateaued 
(Fig.|2]). Here, we report on the highest resolution (1.7 GHz) 
observations which resolve dynamic structures within the flar- 
ing X-ray region, including the discovery of superluminal mo- 
tion in multiple knots. The implications for observations of 
more distant 7-ray emitting relativistic jets are discussed. 

2. DESCRIPTION OF OBSERVATIONS 

Each of our nine VLBA runs consisted of six one-hour inte- 
grations of M87 over an 8 hr period, with calibrator observa- 
tions interleaved. Four adjacent 8 MHz bandwidth channels 
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were used, centered at 1.667 GHz. The full 10 antenna ar- 
ray was used except in Jan. 2005 (missing Hancock antenna). 
May 2005 (missing Saint Croix), and May 2006 (missing Pie 
Town). The data were calibrated using NRAO's VLBA data 
calibration pipeline in AIPS and post-processed with a com- 
bination of DIFMAP and AIPS. The rms noise in the VLBA 
images are within 15% of the theoretically predicted value 
for the full VLBA (0.053 mJy/bm). The images have been re- 
stored with a common beam (Fig. [51) which is within 10% of 
the uniformly weighted beams of all but the May 2005 epoch. 

Over the inner jet, the new sequence of VLBA images shows 
previo usly known radio structure out to "^250 rn as from the 
core ( Biretta& Junoiill995l; iDodson et all l2006h with addi- 
tional transversely resolved emission between ~200-450 mas 
(Fig. [!])■ The central position angle (PA--290°) of the first 
^100 mas portion of the jet aligns roughly with the large- 
scale VLA jet. Knot 'HST-l' is further downstream and is 
offset from the central jet being aligned with the projection 
of the northern edge of the resolved ^200-450 mas jet (i.e., 
PA«294.5°). 

A short 2.3 G Hz VLBA obser vation obtained in Jul. 2004 
by the USNO'' (iFey & Chariot ' 2000) detected only the two 
brightest features in HST-l seen in our early VLBA images al- 
lowing us to trace these structures back by ^1/2 year before 
our program commenced. Further archival VLBI observations 
of va rying quality as far back as in 2000 (e.g., Dodson et ai] 
120061) did not give any significant detections of HST-l. We 
registered the images from the different epochs on the posi- 
tion of the maximum of the radio core ( Fig.[l]l. Based on the 
stability of the structure of the inner jet (iDodson et alJl2006h 
over the 10 total epochs of observations, the core position is 
determined to be aligned to a fraction of the common beam 
over the two year period. 

With the images registered, the upstream (eastern) edge 
which we call 'HST-ld', is the dominant feature in the HST- 
1 complex at the early epochs. The later epoch images re- 
veal the emergence of a radio knot (HST-lc) moving down- 
stream from HST-ld at a rate of 4.48 ±0.42 mas/yr (J3^pp = 
1.14 ± 0.14c) at a PA of 279° and the peak radio surface 

* This research has made use of the United States Naval Observatory Radio 
Reference Frame Image Database (RRFID) 
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Fig. 1 . — Multi-scale images of the M87 jet. Our VLBA images cover the inner ~1" of the jet, as indicated by the corresponding green box in our VLA 3.6cm 
image (from Dec. 2004; upper left inset) of the large-scale jet. The bottom right inset is a zoom-in of the HST-1 complex as imaged with the VLBA (1 pc = 13 
mas); the lowest contour is 0.3 mjy/bm. 



brightness decays by only '^20% over a 1 year period. Some- 
time between Dec. 2005 and Feb. 2006, knot HSTlc evi- 
dently splits into two roughly equally bright features: a faster 
moving component (cl; 4.3±0.7c) and a slower moving trail- 
ing feature (c2; 0.47±0.39c). Between the two years, the lo- 
cation of HST-ld is basically stationary to within r^2 mas (i.e., 
its motion is <0.25c) at 860 mas from the core and is the ap- 
parent point of origin of the superluminal ejections. 

The most distant knot (HST-1 a) is well-isolated from the 
other structures in every observation. We observe it moving 
downstream at 2.49±0.25c at PA = 295°±8° (basically radial) 
until Dec. 2005 when it appears to decelerate to 1.41±0.49c 
at a smaller PA of 289° (although the PA change is apparent 
in Figure in it is not statistically significant). A fainter fea- 
ture (HST-lb), identifiable beginning 2005, trails HST-1 a at 
an identical speed (2.52±0.14c) with a non-radial trajectory 
(a smaller PA of 279°±6°). Feature cl (see above) actually 
ends up in Jul. 2006 where HST-lb was first detected in Jan. 
2005. 

Previous reports of motions up to ~6 c from H ST-1 were 
based on yearly HST optical monitoring jBiretta e t al. 199^ 
from 1994 to 1998 but pertained to emission further down- 
stream from the VLBA structures discussed here. Further- 
more, our motions are apparent on month timescales rather 
than years because of the higher spatial resolution of our 
VLBA data. The extent of the radio emiss ion in HST-1 i s only 
^40 mas, roughly 3 pixels across in the iBiretta et al.l (119991) 
HST FOC images; a Chandra ACIS pixel is ~30x larger. 

3. DISCUSSION AND SUMMARY 

The importance of the superluminal motions reported here 
is twofold. First and foremost, the HST-1 complex is well- 
isolated from the nucleus and the rest of the jet, so the out- 
bursting higher-energy (X-ray, optical) emission isolated by 
Chandra and HST can be uniquely attributed to the region re- 
solved by our VLBA observations (Fig.!!}. While VLBI detec- 
tions of superluminal motions are now com monplace in mor e 
distant relativistic jet sources (e.g.. iKellermann et al.ll2004h . 
the physical link to the higher-energy activity has not been 



possible previously be cause of the lack o f comparably high 
spatial resolution (e.g.. lJorstad et al.ll2001h . Second, contrary 
to conventional wisdom, this 'blazar'-like activity is clearly 
displaced from the central engine by >120 pc (deprojected; 
see § 13.2b . Thus, at least in the case of M87, the observed 
hallmarks of blazar behavior are not directly associated with 
the immediate vic inity of the SMBH where the jet is launched 
dJunor et al.lll999 ). Without the comparably high linear res- 
olution, an analogous sequence of events in a more distant 
source would be associated with the base (i.e., sub-pc to few- 
pc scales) of the jet where the 'non-jetted' contribution from 
the active galactic nucleus (A GN) to the observed X-rays 
(such as from the accretion disk: lMarscher et ani2002h would 
contaminate the lightcurves. 



3.1. The Origin ofVHE j-ray Emission in M87 

At even higher-energies, the H.E.S.S. collaboration recently 
repor ted flaring TeV emission from M87 ( Aharonian et al] 
l2006h . This emission revealed gradual (year timescale) vari- 
ability, with a maximum coinciding with the peak of the 
radio-to-X-ray activity detailed in HST-1 (mid-2005; Fig.|2]l 
suggesting a link between these flares. However, the liST- 
1 knot was dis missed as a possible production site of the 
TeV emission ( Aharonian et al.l [200 6) because of the short 
(days) timescale TeV variability detected on top of the longer 
timescale variability. Such rapid variations were considered 
unlikely for HST-1 since they imply (through the causality ar- 
gument) an emission region size^ Ry^ ^ 0.002 (5 pc. We can 
neither claim nor reject the presence of ~day time-scale vari- 
ability in our optical or X-ray data for HST-1 due to insuffi- 
cient sampling. However, the VLBA data show the compact 
knots in HST-1 to be essentially unresolved with semi-minor 
axes <0.15 pc, and the current optical/X-ray variability data 
constrain /?var ^ 0.022 (5 pc (see below), approaching the size 
limits set by the variability of the TeV emission. 



The Doppler factor is <5 = l/F (1 - Vl-P-^ cos 6»), where T is the bulk 
Lorentz factor, and 6 is the angle to the line of sight. 



Superluminal Motion in HST-1 
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Fig. 2. — Lightcurves of the total TeV emission from M87 (taken from 
|Aharonian et al. 2006) and of the jet knot HST-1 (X-ray and two radio bands). 
The 2 keV and 1 5 GHz data up to Aug. 2005 were previously published in 
[Harris et alj t2006). Our subsequent observations now show the X-ray inten- 
sity of HST-1 decHning steeply, similar to the total TeV emission from M87. 
The 1 .7 GHz VLBA points are integrated from the entire HST- 1 complex. 

Here, we suggest HST-1 as a plausible site for the produc- 
tion of a dominant portion of the detected VHE TeV emis- 
sion due to inevitable inverse-Compton upscattering off am- 
bient photon fields (e.g., starlight) by the electrons producing 
the flaring synchro tron optical-to -X-ray emission. This was 
in fact proposed bv lStawarz et alj (|2006) in discussing earlier 
(2004 and before) TeV detections of M87. The connection 
is strengthened by the fact that the maxima of the total TeV 
flux density from M87 and synchrotron (HST-1 only) flares 
observed in 2005 were coincident, and that their luminosi- 
ties are comparable (Fig. In addition, the 0.4- 10 TeV 
emission from M87 at this time is well described by a single 
power-law with spectral index = 1.2 ±0.15 (Si, oc v'"), 
while the optical-to-X-ray power-law slope of HST-1 during 
the same period is similar: aox = 0.99 ± 0.03. These observa- 
tions are consistent with a common origin for the flaring radio 
th rough 7-ray emission from HST-1 as outlined in more detail 
in lStawarz etall (l2006h . 

The above conclusion is supported by the lack of any 
other plausible production site for the VHE 7-ray emission 
in M87. For example, the innermost (sub-pc scale) jet region 
is char acterized by only small-amplitude optical/X-ray vari- 
ability ([Harris et al. 2006) and only mildly relativistic radio 
features (L v et alji2007i) . As argued also by Aharoni an et alj 
(|2P06), this is inconsistent with models for the gene ration of 
VHE emission in the unresolved core dGeorganopoulos et alj 
l2005h . Instead, lAharonian et alJ (l2006l) c onsidered curva- 
ture radiation of ultra-high energy protons dLevinsonl 120001 : 
iBoldt & Loewensteinll2000 [) accelerated by a strong magnetic 
field in the closest vicinity (^ ^Rg) of the SMBH at the center 
of M87. However, this interpretation is problematic also due 
to the fact that the nearest environments of active SMBHs are 
expected to be opaque to TeV photons due to photon-photon 
(77) annihilation on ambient photon fields such as from the 
accretion disk. 

In M87, the 7-ray photons in the energy range covered 
by H.E.S.S. (e = 1 - 10 TeV) interact mostly with pho- 
tons emitting at ^ 2mlc*/he ~ lO'^"'"* Hz, i.e., infrared 
to optical ones. Quantitatively, the optical depth for the 
77 annihilation process can be approximated by t^^ « 
(1 /3)aT: rno, where ctt is the Thomson cross-section. The 
number density no of the hvo nuclear photons within a spher- 
ical volume around the SMBH of radius r is related to 
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Fig. 3. — Representative sequence of VLBA images of knot HST-1 (see 
Fig. 03 restored with a 8.0 mas X 3.4 mas (0.62 pc X 0.26 pc) beam at PA=3° 
(bottom left). The vertical spacing is proportional to time elapsed and straight 
lines trace the motions of the different features (note the apparent deceleration 
of knot "a"). The first epoch is a shorter 2.3 GHz dataset resulting in the 
excess background; subsequent ones are full track 1.7 GHz observations. 

the monochromatic luminosity Lq by hiyo no = Lq /4tt r^c (the 
main contributor to the photon energy density is judged to 
be from the accretion disk). Hence one can find, r^^ w 

200 (Lo/10'*"ergs-i) (i/q/IO^^Hz)"' {r/Rgy\ The consid- 
ered region would be transparent to the TeV photons only 
if T^~f < 1. However, at r ~ (3 - 10) x Rg we expect 
T^^ ^ 1, since t he obser ved luminosity of the M87 nucleus 
dKharb & Shastri 2004; Whvsong & Antonucc: 
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at i/Q = 

Hz is Lo ^ 3 X lO'*" erg s"' , consistent with ADAF ap- 



proximations of the accretion disk (Di Matteo et al. 2003). 
With the central region of the AGN excluded as a plausible 
site of the TeV emission, the flaring HST-1 knot is the most 
probable point of origin. Therefore, in addition to the flaring 
synchrotron (radio-to-X-ray) emission, the superluminal ra- 
dio structures detailed here can also be associated with VHE 
7-ray activity. 

3.2. HST-1 as a Recollimation Shock 

A natural question to ask is if the position of the HST-1 
knot in the M87 jet is 'special' in some way. In other words, 
what determines the location and flaring behavior of this ex- 
tremely compact and variable feature within the relativistic 
outflow? The kinematics of the superluminal radio compo- 
nents (in particular, the 4.3c motion in cl) constrain the jet at 
the location of HST-1 to be aligned at 6* < 26 ±4°, and re- 
quires r > ^Pl^p+l = 4.4, and S >2. The HST-1 region 
is therefore >120 pc (deprojected assuming 9 = 30°). For 
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Fig. 4. — Positions of radio features in tlie HST-1 complex relative to the 
core (at origin) from our VLBA data (see Fig. [3). The positional errors are 
l/6th of the uniform weighted beam-sizes. Component 'c' emanates from the 
basically stationary feature 'd' (black cross) and splits into cl and c2. All 
components are moving downstream (upper right). 



a SMBH with Mbh « 3 x IO'^Mq (iMacchetto et al. lfT997h . 
this corresponds to the distance of ^10* Schwarzchild radii 
(Rg = GXbh/c^ «4.4 X 10'"* cm « 1.4 x IQ-Vc). At such a 
large distance from the active center, our X-ray observations 
establish flux doubling timescales fvar of 0. 14 yrs (iHarris et alJ 
I2OO6). which constrains /?var < cfyariJ ~ 0.022^pc. 

Recently, Stawarz et al. (2006) showed that the position of 
the HST-1 knot agrees with the expected location of a 're- 
confinement nozzle' formed within the M87 jet due to a con- 
verging shock driven by the interaction of the outflow with 
the interstellar medium (ISM). In this ambient medium, it 
was postulated that the gravitational influence of the central 
SMBH in the inner region of the M87 host galaxy ( Lauer et al.l 
[T9921) results in an increase of the gas pressure, analo- 
gous to the formin g of the observed enhanced stellar cusp 
dYoungetaTII 19781) . For the temperature kTisn, w 0.8 keV of 
hot gaseous ISM (with the number density rii^m ~ 0.17cm~^) 
in the inner (< Ikpc) parts of the galaxy, the gravitational 
capture radius is Ra = G Mbh/c^ ^ lOOpc, wh ere Cs is the 
appropriate sound speed (iDi Matteo et al.ll2003h . It turns out 
that ^ lOOpc is the spatial sc ale of the disk of ionized gas 
observed in M87 by the HST dFord et al.1119941) . It was the 
Keplerian rotation of the inner parts of this disk which en- 
abled a precise estimate of the black hole mass in this system 
( IMacchetto et al. II 19971) . 

Thus, one can surmise that a gravitationally perturbed am- 
bient medium leads naturally to the formation of a jet fea- 
ture like HST-1 at lOOpc from a central SMBH. In this 
scenario, the 'stationary' region HST- Id defines the opening 



of the nozzle. The origins of the different superluminal fea- 
tures can be traced back to this position at ~2004.5 (HST-lc), 
2003.3 (HST-lb), and 2001.9 (HST-la). The latter two could 
have originated more recently as they assume the nominal ve- 
locities of 2.5c; the apparent curve in trajectory and deceler- 
ation of HST-la lead us to believe it could have been moving 
faster previous to our first VLBA imaging observations. 

The formation of a 'reconfinement nozzle' within the 
jet is no t unique for t he hyd rodynamical outflow consid- 
ered by IStawarz et al. (2006). It was previously shown 
dTsinganos & Bo govalov 2002) that interactions between a 
strongly magnetized relativistic outflow with a non-relativistic 
collimating magneto-hydrodynamical (MHD) wind may lead 
as well to the formation of a converging shock within the rela- 
tivistic component. An evaluation of the exact position of this 
shock from the jet base is extremely model-dep endent. Nev- 
erthele ss, for the jet parameters anticipated by iGracia et al.l 
( |2005|) . who successfully explained the observed gradual col- 
limation of the M87 jet between 10^- 10^ Re (ljunor et al.1 
119991) based on the iTsinganos & Bogovalov! (l2002l) MHD 
model, the expected position is near the location of the HST-1 
knot. We speculate that a rapid release of the magnetic en- 
ergy within the formed nozzle via magnetic reconnection may 
play an important role in producing the observed high energy 
flares and ejection of the relativistic radio blobs, somewhat 
analogous to coronal mass ejections from the Sun (Pick et al] 
l2006l) . 

In summary, we have discovered superluminal motions of 
radio features in the M87 jet at a site remote from the central 
SMBH. These features appear to be associated with the re- 
markable flare of 2005 for which the radio, optical, and X-ray 
flux densities peaked at levels of >30-50 times that of a few 
years earlier We have argued that the TeV excess intensity 
of 2005 detected from M87 by H.E.S.S. was a manifestation 
of the same event. Thus, most of the defining characteristics 
of blazars have now been shown to occur at a distance >120 
pc from the SMBH instead of from the immediate environs 
of the central engine, as is commonly believed to be the case. 
As we enter the GLAST era, ensuing studies of high-energy 
flares from blazars should consider such a production site as 
we have resolved in M87. 
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